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UNLIMITED 


SUMMARY 

Since  1975  the  French  have  conducted  underground  nuclear  weapon 
tests  beneath  the  Mururoa  atoll  in  the  Tuamotu  Archipelago  in  the  Southern 
Pacific.  A  least  squares  joint  estimate  of  origin  time  and  epicentre  together  with 
a  least  squares  joint  estimate  of  magnitude  is  presented  for  presumed  explosions 
in  this  area.  These  are  based  on  data  taken  from  bulletins  of  the  International 
Seismological  Centre.  Some  seismograms,  recorded  in  Canada,  are  included  to 
illustrate  the  nature  of  the  P  waves  from  these  presumed  explosions. 


1.  INTRODUCTION 

Basic  source  information  (location,  origin  time,  depth,  yield,  etc) 
about  underground  nuclear  explosions  is  important  to  seismologists  interested  in 
fundamental  properties  of  the  structure  of  the  earth,  as  well  as  those  interested 
in  discrimination  between  earthquake  and  explosion  from  their  seismic  signals. 
Numerous  scientists  have  appealed  for  the  release  of  epicentral  details  of 
explosions  to  aid  research  programmes  (Bullen,  1958,  Griggs  and  Press,  1961, 
Teller,  1963).  In  response.  Springer  and  Kinnaman  (1971,  1975)  published  the  basic 
epicentre  details  for  all  underground  nuclear  explosions,  US  and  British, 
detonated  in  the  USA  from  1961  to  1973.  Numerous  yield  estimates  were  also 
included.  The  origin  time  and  precise  epicentres  of  French  underground  nuclear 
explosions  in  the  Sahara  between  1961  and  1966  have  been  published  by  Duclaux 
and  Michaud  (1970).  However  no  comparable  data  are  available  for  underground 
explosions  in  the  USSR  or  the  French  Pacific  Test  Site  at  the  Tuamotu 
Archipelago. 

Several  international  data  centres  collect  seismic  wave  arrival  times 
from  all  over  the  world  and  compute  estitnates  of  the  origin  time,  epicentre, 
depth  and  size  for  seismic  disturbances.  Bulletins  containing  these  data  are 
published  by  the  US  National  Earthquake  Information  Centre  (NEIC)  in  Colorado, 
USA,  and  the  International  Seismological  Centre  (ISC)  in  Newbury,  UK.  A  similar 
service  is  provided  by  the  Institute  of  Physics  of  the  Earth  in  Moscow,  but  the 
Soviet  bulletin  does  not  report  data  on  Soviet  nuclear  explosions.  The  French 
produce  a  bulletin  for  their  metropolitan  and  Polynesian  seismological  networks 
but  the  bulletin  contains  no  information  on  French  explosions. 

To  provide  seismologists  with  improved  estimates  of  some  of  the 
basic  source  information  about  French  presumed  underground  nuclear  explosions, 
seismic  wave  arrival  time  and  amplitude  data  collected  by  the  ISC  have  been 
analysed  using  a  joint  epicentre  technique  (JED)  to  relocate  the  epicentres,  and  a 
least  squares  analysis  (LSMF)  of  amplitude  data  to  provide  consistent  estimates 
of  the  seismic  magnitude.  A  similar  study  has  already  been  made  for  explosions 
at  the  Soviet  test  site  near  Semipalatinsk,  East  Kazakh  and  has  s^en  published 
by  Marshall  et  al  (198^*). 


It  is  planned  to  carry  out  a  similar  study  of  other  Soviet  test  sites 
such  as  Novaya  Zemlya,  North  Caspian  Sea  region  and  the  Nevada  Test  Site 
(NTS)  USA. 

2.  JOINT  EPICENTRE  RELOCATIONS 

The  Joint  Epicentre  Determination  (JED)  method  described  by 
Douglas  (1967)  was  employed  to  relocate  the  events  using  P  &  PKP  arrival  time 
data  taken  from  International  Seismological  Centre  (ISC)  bulletins.  To  fix  the 
overall  location  of  the  group,  one  of  the  epicentres  must  be  restrained  to  a 
predetermined  value.  The  restrained  epicentre  chosen  was  that  for  the  event  on 
25  July  1979,  the  largest  and  most  widely  recorded.  No  true  epicentre  has  been 
published  and  so  the  location  must  be  fixed  using  other  evidence.  The  strategy 
used  here  was  to  shift  the  restrained  epicentre  until  the  overall  pattern  fitted 
centrally  over  the  island  (figure  I).  Taking  the  median  latitude  and  longitude  as 
representative  of  the  group,  only  when  the  restrained  epicentre  is  located  near 
the  southern  coast  does  this  median  lie  on  the  island  and  its  central  lagoon.  The 
adopted  epicentre  (2I.88S,  138.94W)  gives  the  minimum  deviation  of  the  median 
location  from  the  lagoon  centre  (taken  as  21.83S,  138.91W).  Confidence  in  the 
restrained  location  is  gained  by  the  fact  that  the  ISC  location  (21.86S,  1  39.0W)  is 
approximately  6  km  WNW  of  the  chosen  position,  a  bias  similar  to  that  expected 
when  station  travel  time  corrections  are  not  used  (Lilwall  and  Underwood,  1970). 

Arrival  time  readings  were  weighted  to  remove  gross  errors  and  allow 
for  variation  in  the  quality  of  arrival  time  measurement  between  stations.  The 
effects  of  gross  errors  were  reduced  using  the  method  of  uniform  reduction 
(Jeffreys,  1961).  Where  a  station  reported  sufficient  events  (here  set  at  10)  the 
standard  deviation  of  the  residuals  was  calculated  and  the  station's  arrival  time 
readings  weighted  to  account  for  variation  in  this.  This  technique  also  permitted 
the  incorporation  of  a  large  body  of  PKP  data  which  would  all  normally  be  given 
zero  weight. 

Table  1  gives  the  relocated  epicentres,  origin  times  and  the  95% 
confidence  limits.  All  depths  were  constrained  to  a  surface  focus  source. 
Figure  1  shows  the  epicentres  plotted  on  a  map  of  the  atoll.  Although  the 
standard  deviation  of  the  residuals  (0.28s)  is  small  the  poor  azimuthal  distribution 
of  stations  results  in  relatively  large  limits,  the  majority  (22  out  of  28)  however 
are  less  than  10  km.  All  but  three  (events  11,  22,  24)  have  95%  limits  which 
intersect  the  island. 

3.  DETERMINATION  OF  MAGNITUDE 

The  size  of  a  seismic  source  is  measured  by  its  magnitude.  For  short 
period  (SP)  seismic  P  wave  data  the  Gutenberg  and  Richter  definition  is  used:- 

=  logjjA/T  +•  B(A)  S,  ...  (1) 

where  A  is  the  amplitude  of  the  P  wave  in  nm,  T  its  predominant  period  in 
seconds,  B(A)  a  distance  normalising  term  and  S  a  station  correction. 


Consider  n  explosions  recorded  at  some  or  all  of  q  stations.  Then  if  m 
is  the  magnitude  of  the  ith  explosion  recorded  at  station  j,  we  can  write 


=  b .  +  s  . 
1  j 


where  b.  depends  on  the  seismic  size  of  the  explosion,  s.  is  a  station  term  and  e .. 
is  an  error  term.  Least  squares  can  be  used  to  estimate*  b;  and  S;  if  it  is  assumed 


By  making  the  further  assumption  that  the  errors  e  ..  are  normally  distributed, 
confidence  limits  can  be  determined  for  b.  and  s..  Tms  is  the  procedure  used  to 
obtain  the  least  squares  estimate  of  the  magnitude  of  explosions  (Douglas,  1966). 

The  basic  input  data  are  taken  from  the  ISC  bulletins,  either  as  m.  or 
in  the  form  of  Jog^^A/T  from  which  fT>.  is  calculated.  The  stations  used  are 
located  in  the  distance  range  60  to  90  ."The  mean  magnitude,  (a  least  squares 
estimate)  of  the  explosions  is  given  in  table  1.  Station  corrections  are  given  in 
table  2  together  with  the  95%  confidence  limits  and  the  number  of  observations. 
The  statistics  associated  with  the  least  squares  analysis  are  given  in  table  3. 

U.  SEISMOGRAMS 

The  short  period  seismological  array  station  located  at  Yellowknife, 
Canada,  appears  to  be  well  placed  to  detect  explosions  at  Mururoa.  Several 
examples  of  the  beamed-array  P  wave  are  illustrated  in  figure  2.  It  would  appear 
that  the  seismograms  fall  into  one  of  two  categories.  For  example,  explosion 
number  19  has  no  large  second  trough  whereas  explosion  number  24  has  a  clear, 
large  second  trough;  the  presence  or  absence  of  the  second  large  amplitude 
trough  defines  the  two  categories.  The  presence  of  the  second  trough  may  be  a 
"depth"  effect  caused  by  the  interaction  of  the  direct  P  with  the  free-surface 
reflection  (pP)  since  most  of  the  explosions  which  exhibit  this  feature  are  of 
relatively  large  magnitude  compared  to  the  lower  magnitude  explosions  for  which 
the  second  trough  is  absent. 
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TABLE  ?(3j 


Station  Terrs  for  Explosions  at  Mururoa 
ITE21ATI0N  0 


STATION 

COMPUTED* 

VALUE 

9?  PERCENT 

coih7ide:;ce  limits 

NUMBER 

kUi 

-O.P^L 

V-  c.ll? 

16 

ARE 

-O.P17 

V-  0.280 

2 

ASP 

-0. 

V-  0.229 

3 

3CT 

0.09? 

■>•/-  0.288 

2 

BDW 

-C. JA? 

V-  0.158 

7 

B7D 

-0.03? 

V-  C.20C 

9 

BKS 

0.?lU 

V-  0.129 

10 

BU 

-0.907 

+/-  0.229 

3 

BMN 

-0.190 

V-  0.107 

18 

BNH 

-0.097 

V-  0.280 

2 

COL 

0.090 

■*/-  0.121 

12 

CTA 

0.197 

-/-  C.20c 

9 

DUG 

-O.l'fC 

-/-  0.?^9 

T 

EDM 

0.69^ 

■*■/  -  1^- 

8 

EMM 

0.33? 

T  2  ’  A’ 

5 

EUR 

-C.9c- 

+/_  C.I27 

11 

FBA 

O.OEl 

*•/-  0,1='^ 

7 

FFC 

-c.-oc 

V-  C.2-9 

■* 

rvM 

0.63c 

V-  C-?’3 

7 

GIL 

0.CC.1 

V-  C,215 

9 

GLS 

-O.C':-' 

V-  C-285 

2 

hjCirf 

-0.951 

v-  c.i’.: 

9 

KDM 

0.3?r 

-/-  0.23'' 

7 

HNH 

0.679 

+/-  0.282 

2 

JCT 

-0.099 

V-  0.190 

3 

1I'3 

-O.C99 

V-  0.183 

c 

J 

LON 

-0.086 

V-  0.200 

9 

LP5 

-C.?C8 

V-  o.if  ’ 

6 

LPS 

0.06’ 

V-  0.207 

9 

MAW 

C.?5? 

V-  0.199 

9 

MSO 

0,169 

V-  0.105 

le 

NEW 

O.C’i* 

V-  C.2C1 

*9 

NNA 

-O.A*-? 

V-  C-78c 

2 

PI-IG 

-0. 029 

‘/-  C.l'’8 

c 

P>K 

O.c^C 

V-  0.097 

19 

PNT 

0.  Cr  8 

V-  C.l:c 

5 

SDV 

C.3-8 

+/ —  0. 200 

9 

SES 

0,3=? 

♦/-  C.I99 

9 

SPA 

-0.58-; 

-/-  C.?OC' 

H 

STV. 

C.^0^5 

V-  0.2 ’1 

T 

TOC 

C.-3’ 

-/-  0.229 

7 

TCV 

-0.C?5 

>/-  0.2^0 

TUC 

-0.030 

7 

TUL 

0.1  ?6 

+/-  0.119 

12 

UAV 

C.?45 

-  G .  V  30 

7 

UCT 

C.  ?”C 

-/-  0.281 

? 

WRA 

-0.-8? 

T.  -  0.10*t 

19 

Y};c 

0.?5- 

V-  o.iit 

14 

*  To  u?e  as 

staticr.  correrticr.s  the  sien  of 

the  station 

terr.  must 

be  reverse! 

(see  e-nuation  (/)). 
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TABLE  3 


Statistical  and  Other  Variables  Used 


Iteration  0 


Number  of  Readings 
Number  of  Unknowns 

Sum  of  Squares  of  Residuals 
Mean  Square  of  Residuals 
Number  of  Degrees  of  Freedom 
Students  T 
Sum  of  Residuals 

Sum  of  Squares  Due  to  Stations 
Mean  Square  Due  to  Stations 
Number  of  Effective  Stations 

Sum  of  Squares  Due  to  Events 
Mean  Square  Due  to  Events 
Number  of  Effective  Events 


French  Test  Site 


=  296 
=  76 

=  8.923898 
=  0.040563 
=  220 
=  1.97 
=  0.36E-12 

=  23.121092 
=  0.491938 
=  47 

=  33.634553 
=  1.201  234 
=  28 


Mean  Station-Event  Effect 
Sum  of  Dummy  Variables 
Squares  Due  to  Dummy  Variables 


5.112  ±  0.04483 

0.32E-14 

0.84E-45 


TABLE  3  (continued) 

Iteration  3 

Number  of  Readings 

s 

296 

Number  of  Unknowns 

s 

76 

Sum  of  Squares  of  Residuals 

3 

2.392113 

Mean  Square  of  Residuals 

s 

0.010873 

Number  of  Degrees  of  Freedom 

= 

220 

Students  T 

= 

1.97 

Sum  of  Residuals 

= 

-0.80E-03 

Sum  of  Squares  due  to  Stations 

= 

19.978145 

Mean  Square  due  to  Stations 

0.425067 

Number  of  Effective  Stations 

= 

47 

Sum  of  Squares  due  to  Events 

s 

33.954501 

Mean  Square  due  to  Events 

ss 

1.212661 

Number  of  Effective  Events 

= 

28 

Mean  Station-Event  Effect 

5.100  ±  0.02369 

Sura  of  Dummy  Variables 

= 

-0.32E-04 

Squares  due  to  Dummy  Variables 

s 

0.22E-15 

FIGURE  1.  PLOT  OF  COMPUTED  EPICENTRES  ON  A  MAP  OF  THE  MURUROA 
ATOLL.  (NOTE  THAT  EVENTS  NUMI^EREO  10,  II  AND  22  IN 
TABLE  1  ARE  LOCATED  OUTSIDE  AREA  OF  MAP). 
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FIGURE  2.  EXAMPLES  OF  THE  P-WAVEFORMS  RECORDED  AT  YELLOW¬ 
KNIFE,  CANADA  FROM  UNDERGROUND  EXPLOSIONS  AT 
MURUROA. 
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